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Abstract — We develop qubit Hamiltonian single parameter 
estimation techniques using a Bayesian approach. The algorithms 
considered are restricted to projective measurements in a fixed 
basis, and are derived under the assumption that the qubit mea- 
surement is much slower than the characteristic qubit evolution. 
We optimize a non-adaptive algorithm using particle swarm 
optimization (PSO) and compare with a previously-developed 
locally-optimal scheme. 

I. Introduction 

The field of quantum computation is one of the fastest 
developing areas in modern physics. In order to perform 
quantum computation one needs to have a full and precise 
information about the parameters of quantum bit evolution in 
the set-up used. This makes quantum Hamiltonian parameter 
estimation one of the crucial problems in this area. 

The most common approach involves repeatedly initiating 
the system in some state and then performing a measurement 
after some time when the system is allowed to evolve under 
its Hamiltonian (TJ. The choice of input states, evolution 
times, measurement bases and method of information treat- 
ment defines the estimation algorithm. Various approaches 
to the problem were suggested, including quantum process 
tomography [2|, [3 | and Fourier analysis [4|. Quantum process 
tomography requires a full control over state preparation and 
measurement which could be a hard task. In addition, often 
the preparation is done by using strong measurement in a 
certain basis which one can change only using evolution of the 
quantum bit (qubit) itself which is impossible without already 
knowing the Hamiltonian parameters. Thus finding algorithms 
which are able to work using a fixed basis is required. Fourier 
analysis method is inspired by classical signal processing. 
It involves repeated measurements after longer and longer 
evolution times and signal frequency reconstruction using 
inverse Fourier transform. However, it requires a large number 
of measurements to achieve moderate precision. I.e., in |0] it 
took 10 6 measurements to achieve 10~ 3 joint variance in a 
two-parameter qubit Hamiltonian estimation. 

In order to completely use the information of each measure- 
ment it is necessary to use a Bayesian approach. It consists of 
representing the knowledge of an unknown parameter in terms 
of probability distribution and updating it according to the 
sequence of measurement outcomes. This approach is closely 
related to phase estimation techniques, which however also 
employ arbitrary states preparation JT], |]5]. 



The particular problem considered in this paper is inspired 
by the developments in 2-electron double quantum dot (DQD) 
systems which can be successfully used for creating quantum 
bits and performing gate operations, in particular GaAs double 
quantum dots J6), 0. In such systems qubit evolution is 
driven by a difference in z-component of magnetic field AB 
between the parts of the DQD. This field is due to the hyperfine 
interaction of the electrons' spins with atoms nearby lattice and 
so it is not stable and can't be controlled. Precise and quick 
measurement of this field is therefore crucial for quantum 
processing. One of the important limitations imposed by this 
system is a slow and complicated measurement process which 
requires considerably more time than the evolution (~ 6/xs 
vs several nanoseconds). In order to optimize the estimation 
procedure we'll quantify the resources which are needed to be 
minimized in terms of number of measurements (not the total 
evolution time as in previous research available, i.e. ||8]). 

The optimization process for such algorithm is hard as 
in general it requires sampling the probability distribution 
and browsing through an enormous amount of algorithm 
variants. The dependence of the final estimate variance (to 
be minimized) on the algorithm is a complicated function 
with many local minima and unpredictable behaviour. Particle 
swarm optimization (PSO) approach has proved to be useful 
in such situations when little is known about the function and 
global minimum is sought iflOl . ifTTl . It was also recently used 
in quantum phase algorithm optimization lfl2l . Ifl3l . In this 
paper we present a locally-optimal approach derived in |]9] 
and demonstrate the improvement which can be achieved by 
using particle swarm optimization technique. 

II. Hamiltonian parameter estimation 

A. Problem statement 

We consider the problem of single unknown parameter 
estimation of a qubit Hamiltonian of the form 



H 
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(1) 



where a z is a z-Pauli matrix. 

The system is initiated in {| — ) , |+)} basis (also referred to 
as "computational basis") by applying a projective measure- 
ment along z-axis. The time evolution of the state initiated as 
|— ) under the Hamiltonian (HJ can be represented as 



\ip{t)) = cos (wi/2) sin (wt/2) |+) 



(2) 



The measurement in the computational basis is then performed 
giving either |— ) or |+) outcome. 

We assume that the unknown parameter lu can be treated as 
a random variable uniformly distributed in (0,luq). The upper 
bound is defined by Nyquist-Shannon theorem as luq = n/ At 
where At is the minimal accessible time interval. Using time 
intervals quantized as t = mAt, m e N it is easy to see 
that the evolution produces the following probabilities of the 
measurement outcomes: 
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Assuming that each measurement takes the same amount 
of time, the problem would become to optimally select the 
sequence of evolution times m^At. 

B. Bayesian approach 

Treating the information about w as a probability distribu- 
tion we now can use Bayesian updating to find the estimate. 
Taking the flat prior, distribution after some number of steps 
will merely be a normalized product of probability functions 
(01. The variance of this distribution can thus characterize the 
precision in parameter estimate. 

The choice of discrete evolution time selection as discussed 
above is essential for brute force algorithm optimization. In 
such case we can avoid sampling the distribution by repre- 
senting it as a Fourier series. It can be done for any algorithm 
as p(±\ui) are both harmonic functions. After k measurements 
the distribution reads 

I M 

Pk{uj\r k ■ ■ ■ n) = -Cfc(O) + J2 c k(q) cos(gvrw/w ), (4) 

^ q=l 



and 



are the measurement outcomes. 



where M ~ m j 

i=i 

The distribution is normalized by dividing by ^Lu Ck(0)- The 

lu 2 ) — (lu) 2 can then be 



variance of this distribution V 
found as 
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The expected variance E[V] is now the function to minimize 
and it depends only on the selection of the set of evolution 
times. Also it is clear that E[V) is constant under the permu- 
tations of elements of {rrii} set. 

III. Algorithm optimization 

A. Locally-optimal approach 

The straight way to optimize the algorithm is to use local 
optimization |9), lfl4l . This approach involves fixing first k 
evolution times and optimizing the expected variance after the 
next measurement over the choice of irik+i- Repeating this 
procedure for k = 1,...,N we create a locally-optimal non- 
adaptive (LONA) algorithm. First 10 steps of LONA are shown 
in Table [TTJ The performance of LONA algorithm is shown on 
Fig. [T] as a solid line. 



TABLE I: LONA and optimized algorithm for 5 measure- 
ments. 



Step 


LONA 


PSO 


1 


1 


1.060 


2 


1 


1.082 


3 


1 


1.419 


4 


2 


2.138 


5 


3 


2.870 



TABLE II: LONA and optimized algorithm for 10 measure- 
ments. 



Step 


LONA 


PSO 


1 




1.071 


2 




1.107 


3 




1.161 


4 




1.180 


5 




1.200 


6 


2 


2.041 


7 


2 


2.152 


8 


3 


3.070 


9 


4 


3.970 


10 


6 


4.906 



As it was argued in (9), this approach demonstrates a 
significant improvement over a naive approach with fixed 
evolution time rrik = 1 and Fourier techniques. However, 
numerical simulations show that this approach is not globally- 
optimal being slightly outperformed by mk = k algorithm 
in a short range around 60th measurement. In addition, as it 
was mentioned before, the general algorithm optimized over 
TOfe £ K would be very complex computationally. Firstly, 
it will require distribution sampling, that is instead of a 
set of Fourier coefficients (which perfectly reconstruct the 
function), the distribution should be stored point-by-point with 
limited precision. Secondly, the brute force browsing through 
evolution times at each step and performing the averaging 
will make it unrealistic even for few steps as the number of 
operations grows as a very steep exponent. Globally-optimal 
optimization and selecting all the possible sets of {m^} that 
way becomes completely unrealistic. 

B. Particle swarm optimization 

The globally-optimal non-adaptive technique can be found 
with high precision employing PSO approach. It enables to 
perform the search in all the continuous space of the set 
of evolution times {ti}. We are sampling the lu probability 
distribution using 10 4 points which gives < 10~ 6 precision in 
E[V] for up to N = 10. 

Since the evolution times are permutable, we need to avoid 
covering the identical algorithms. In order to do this the 
function is minimized over the variables dt±, dt2, ...,dtN, 
where dt\ = t\ is the first evolution time and dt k = t k — tk-i- 
The PSO 'particles' thus have iV-dimensional coordinates and 
the 'goodness' of their position is defined as E[V] for these 
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Measurement step 

Fig. 1: Performance of LONA algorithm (solid line) and 
improved algorithms found using PSO (crosses). 




PSO step 

Fig. 2: Joint variance convergence with the number of PSO 
updates for 5-measurement estimation algorithm. 



evolution times. If dtk < we manually set the value of E[V] 
to be some large number (instead of calculating the expected 
variance), creating a barrier for the particles so that they cover 
only unique algorithms. 

We use a swarm of 16 particles, randomly initiated with 
dti £ (0,1.5) and dti £ (0,2), which is a reasonable 
range implied from the discrete LONA algorithm. The initial 
velocities are arbitrarily set to (—1,1) range and their absolute 
values are limited from above by 2. 

We employ the velocity updating rule with a constriction 
factor \ as follows 



where cp — ci + c 2 , as was suggested in |[T5l . Numeric 
simulations showed that this technique performs much better 
than the standard approach with \ an d c i = c 2 = 2, having 
a relatively quick and stable convergence and no 'swarm 
explosion' occurrences. 

The optimization was done for 5 and 10 measurement steps 
using the swarm of 16 particles. For 5 steps we initiated the 
coordinates as being random ±0.5 within LONA algorithm 
and for 10 steps ±0.1. A typical plot for the convergence 
of joint variance of the swarm is shown on Fig. [2] for 5- 
measurement algorithm. The best algorithms found are shown 
in Tables Q] and [TT] and their performance is shown as crosses 
on Fig. [T] 

IV. Conclusion 

We demonstrated how the performance of qubit Hamil- 
tonian estimation algorithm can be improved using PSO. 
This approach requires several orders less operations in order 
to achieve similar precision in evolution times defining the 
algorithm. Thus, we demonstrated how PSO can be applied 
to complex quantum algorithms optimization to achieve no- 
ticeable results and that it can be useful in such applications. 
The further advancements would include creating a realistic 
estimation process model for DQD system in order to optimize 
the procedure for a real experimental set-up. 
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'Ji+1 



X(vi +ric 1 (x l -x) + r 2 c 2 (x g - x)) 



(6) 



where xi and x g are local and global best coordinates so far 
and r%,2 £ (0, 1) are random numbers generated independently 
for each particle every swarm updating step. The parameters 
used are c\ = c 2 = 2.05 and 

X 2 . « 0.729, (7) 



